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It is known that a nucleus with charge Ze, where Z > 170 creates electron-positron pairs from the 
vacnum. Electrons collapse onto the nuclens resulting in a net charge Zn < Z while the positrons are 
emitted. This effect is dne to the relativistic dispersion law. The same reason leads to the collapse 
of electrons to the donor cluster with a large charge number Z in narrow-band gap semiconductors, 
Weyl semimetals and graphene. In this paper, a similar effect of electron collapse and charge 
renormalization is found for a donor cluster in SrTiOa (STO), but with a different origin. At low 
temperatures, STO has an enormously large dielectric constant and the nonlinear dielectric response 
becomes dominant when the electric field is still small. This leads to the collapse of electrons into 
a charged spherical donor cluster with radius R when its total charge nnmber Z exceeds a critical 
value Zc — R/a where a is the lattice constant. The net charge ZnC grows with Z until Z exceeds 
Z* ~ (R/a)®'^^. After this point, the charge of the compact core Zn remains ~ Z*, while the 
rest Z* electrons form a sparse Thomas-Fermi electron atmosphere aronnd it. We show that the 
thermal ionization of such two-scale atoms easily strips the outer atmosphere while the inner core 
remains preserved. We extend onr resnlts to the case of long cylindrical clnsters. We discnss how 
our predictions can be tested by measuring conductivity of chain of discs of charge on the STO 
surface. 


I. INTRODUCTION 

Recently, studies of ABO 3 perovskite crystals have 
been a subject of interest due to their intriguing mag¬ 
netic, superconducting, and multiferroic propertied and 
subsequent signihcant technological applications. Am ong 
them, SrTiOa (STO) has attracted special attentioiPd 
STO is a semiconductor with a band gap of ~ 3.2 eV 
and a large dielectric constant k = 2 x 10“^ at liquid he¬ 
lium temperature. Like the conventional semiconduct ors, 
STO can be used as the basis for a number of deviced^. 

Many of the devices are realized by doping the bulk 
STO such as the reduction of STO through generating 
oxygen vacancies at high temperatures. The vacancies 
either form along a network of extended defectd or as¬ 
semble together to lower the system’s energjl^, prob¬ 
ably producing large positively charged donor clusters. 
Another way to more controllably create such a cluster 
is to “draw” a disc of charge by the atomic force micro¬ 
scope (AFM) tip on the surface of LAO/STO structure 
with the subcritical thickness for LaAlOa (LAO)P. The 
potential caused by such a positive disc in the bulk STO 
is similar to that of a charged sphere. 

Let’s consider a spherical donor cluster with radius R 
and charge Ze. We assume that the cluster is located 
on the background of uniformly n-type doped STO in 
which the Fermi level is very close to the bottom of the 
conduction band. There are Z electrons located at dis¬ 
tances from KbjZ to the Bohr radius nh from the clus¬ 
ter, which form a Thomas-Fermi “atom”Ii^ with it. Here 
b = f? jm*^, m* « l.Srrie is the effective electron mass 
in sTcmi 

with TOe being the electron mass. Since k is 
large, the electrons are far away from the cluster and the 
whole “atom” is very big. As Z increases, the electron gas 
swells inward to hold more electrons. However, we find 


that when Z goes beyond a certain value Zc {nh/Z is still 
much larger than R at this moment), the physical picture 
is qualitatively altered. Surrounding electrons start to 
collapse into the cluster and the net cluster charge gets 
renormalized from Ze to Z„e with Zn Z ai very large 
Z. 

The effect of charge renormalization is not nevn^ESl. 
For a highly charged nucleus with charge Ze, the vacuum 
is predicted to be unstable against creation of electron- 
positron pairs, resulting in a collapse of electrons onto 
the nucleus with positrons emittecP^. This instability 
happens when Z > Zc with Zc — 170 > l/o:, where 
a = e^/he ~ 1/137 is the hne structure constant. When 
Z exceeds ~ ~ 137^/^, the net charge satu¬ 
rates at In the condensed matter setting, there are 

similar phenomena in narrow-band gap semiconductors 
and Weyl semimetalJi^as well as graphen^^l these 
cases, the collapse happens because the energy dispersion 
of electrons is relativistic in the Coulomb field of a com¬ 
pact donor cluster playing the role of a nucleus. In our 
work, however, the collapse originates from the strong 
nonlinearity of dielectric constant in STO at small dis¬ 
tances from the cluster. In the case of a spherical donor 
cluster, this nonlinearity leads to the change of the at¬ 
tractive potential near the cluster from being oc l/r to 
oc l/r^, resulting in the collapse of electrons to the clus¬ 
ter. 

The phenomena of electron collapse and charge renor¬ 
malization in both heavy nuclei and our work are pre¬ 
sented in Fig. In our case, the first electron collapses 
aX Z Zc — R/a, where a is the lattice constant, and 
at Z Z* cx (R/a)®/"^, the net charge of nucleus Z„e 
saturates as Zn — Z*. 

In the remainder of this paper, we use the Thomas- 
Fermi approximation to show how the electron gas col- 
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FIG. 1. (Color online) The number of collapsed electrons 
S and the renormalized net charge Z„e as a function of the 
original charge Ze. S is shown by the thin solid line (red), 
Zn is denoted by the thick solid line (blue), and the dashed 
line (black) is a guide-to-eye where Z„ = Z. Zc denotes the 
critical value where electrons begin to collapse and Z* is the 
saturation point where Z„ stops growing, (a) Collapse of elec¬ 
trons and charge renormalization for highly charged nuclei. 
S' oc Z® at Zc ^ Z <C Z"^. (b) Collapse of electrons and 
charge renormalization for spherical donor clusters in STO. 
ScxZ^^'^ at Zc < Z < Z*. 


lapses into the cluster at Z ^ Z^. and find the corre¬ 
sponding electron density. In Sec. II, we demonstrate 
that when Zc Z <C Z*, the charge renormalization is 
relatively weak and the final charge number Z„ is just a 
little bit below Z. On the other hand, when Z ^ Z*, 
the renormalization is very strong and Z„ is maintained 
at the level of Z*. The inner core of the cluster atom 
with charge Z* is surrounded by a sparse Thomas-Fermi 
atmosphere of Z* electrons in which the large linear di¬ 
electric constant plays the main role. In Sec. Ill, we 
generalize our studies to the cylindrical donor clusters 
and introduce the notion of maximum linear charge den¬ 
sity r]* similar to Z* which is got when the bare charge 
density of the cluster 77 is large. In Sec. IV we consider 
the thermal ionization of donor cluster atoms and ar¬ 
rive at the conclusion that for both spheres and cylinders 
the external electron atmosphere is easily ionized while 
the inner core with charge Z* or linear charge density 
rj* is robust against the ionization. We suggest experi¬ 
mental verification of our theory for a periodic chain of 
charge discs created by methods of Ref. |S1 If neighboring 
disc atoms overlap via the outer electron atmosphere (the 
“bridges”), raising temperature may ionize these bridges 
and sharply reduce the chain conductance. We conclude 
our work in Sec. V. 


II. SPHERICAL DONOR CLUSTERS 


A. Nonlinear Dielectric Response 


It is well known that STO is a quantum paraelectric 
where the onset of ferroelectric order is suppressed by 
quantum fluctuation^^. In such ferroelectric-like mate¬ 
rials, the free energy density F can be expressed by the 
Landau-Ginzburg theory as a power-series expansion of 
the polarization 

F = Fo + -b - EP 

“ 2 4P2 

(in Gaussian units), where Fg stands for the free energy 
density at P = 0 , P is the electric held, r = 47 r/(K— 1 ) ~ 
47 r/K is the inverse susceptibility, k ;i> 1 is the dielectric 
constant, A « O.fP^, Pq = e/a^, a = 3.9 A is the lattice 
constant. We neglect the gradient terms of polarization 
since they play a minor role in the nonlinear regime. We 
assume that the dielectric response is isotropic. This is 
justihed by the small anisotropy in the nonlinear response 
of STCC^. The crystal polarization P is determined by 
minimizing the free energy density P in the presence of 
the electric held E, i.e., SF/SP = 0. This gives 

47rP AP3 

For very big k, the nonlinear term in this expression is 
likely to dominate over the linear one even at relatively 
small P as long as P » in/ kAPq. At distances not 
too far from the cluster, this relationship can easily be 
satished and the dielectric response becomes nonlineai® 


E = 


AP^ 

p2 ■ 


( 1 ) 


Since Pq = e/a^ is really big, we can expect that P < 
in/AP q always holds, which gives E < inP. Then, 


D = E + AttP « 47rP. (2) 


According to the Gauss’s law in dielectric media, we know 
V • P = 47 rp where p is the free charge density. Together 
with Eqs. ([^ and (|^, we have 

V • (V0)i/3 = - ' p (3) 

where (p is the electric potential. 

Below, we show that this nonlinear dielectric response 
leads to the collapse of electrons into the spherical donor 
cluster inside STO when the cluster charge is large 
enough. 


B. Renormalization of Charge 

Gonsider a large spherical donor cluster of the radius 
R and the total positive charge Ze such that a <C P < 
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KbjZ. If the dielectric response is linear, the electrons are 
mainly located at distances between ri = KbjZ and = 
nb from the clustei^i^ where b = lT?/m*e^ « 0.29 ^4^^. 
For a very large k, these radii are huge (r^ = 700 nm 
in STO at liquid helium temperature where k = 20000) 
and the electrons are far away from the cluster. However, 
at small distances, the dielectric response is nonlinear 
and changes the potential form. If the potential energy 
outweighs the kinetic energy, electrons are attracted to 
the cluster and renormalize the net charge. To see when 
this will happen, we look at the specific form of electric 
potential in this situation. We can calculate the potential 
from Eq. ([^. But due to the simple charge distribution 
here, we can get it in an easier way. At r > i? where r is 
the distance from the cluster center, the sphere looks like 
a point charge and D{r) = Zejr^. Using this together 
with Eqs. Q and ([^, one can calculate the electric field 
and get the electric potential as: 

, A fZe\^ 1 

= R<r«n (4) 

with (j){r = oo) defined as zero. Inside the cluster at 
r < R, since the charge is uniformly distributed over the 
sphere, the total positive charge enclosed in the sphere 
of radius r is equal to Zer^/B?, so D[r) = Zer/R?. One 
then gets the corresponding potential 




A fZeY ( 9 1 


4i?V 


0 < r < R 


( 5 ) 


using the boundary condition ^(r = R~) = 4>{r = i?+). 
A schematic graph of the potential energy U (r) = —e(j)(r) 
is shown in Fig. [^by the thick solid line (blue). 

The Hamiltonian for a single electron is H = /2m* — 

e(f>(r), where p is the momentum of the electron and m* is 
the effective electron mass in STCff^. If we approximately 
set p ~ h/2r, we get a positive total energy of the electron 
everywhere when Z is very small. This means there are 
no bound states of electron in the cluster. However, when 
Z is big enough so that Z > Z^, the electron can have 
negative total energy at r < i? and will collapse into the 
cluster. Using Eqs. Q and ([^, we find 

Z, =» ...« »1. (6) 

a a 

As Z continues increasing, more and more electrons get 
inside the cluster filling it from the center where the po¬ 
tential energy is lowest (see Fig. [^. The single-electron 
picture no longer applies. Instead, we use the Thomas- 
Fermi approximatiorP^ with the electrochemical poten¬ 
tial /i = 0, which gives: 


n{r) 


Cl 

63 




3/2 


( 7 ) 


where n(r) is the electron density at radius r, ci = 
23 / 2/3712 « 0 . 1 . 


We assume that the bulk STO is a heavily doped semi¬ 
conductor in which the Fermi level lies in the conduction 
band. On the other hand, due to the relatively high ef¬ 
fective electron mass, the Fermi energy is much smaller 
than the depth of the potential well shown in Fig. [^and 
is thus ignored. Since at r = oo the electric potential 
(j){r) is defined as 0, we then have the electrochemical 
potential /t ~ 0. 

When the number of collapsed electrons S is small, 
their influence on the electric potential is weak. One can 
still use Eqs. 0 and 0 for (j){r) and get the correspond¬ 
ing expression of n(r). At r > i?, since (p{r) is oc I/r^, 
we get n{r) oc I/r^®/^. In this way, we calculate S as 

roo f 7 \ ^/2 

S = J n{r)A'Kr‘^dr = Q.bZ ^ (8) 


where 




■47r(6/Aa)3/3i? 

a 



( 9 ) 


The net charge number of the cluster is 


Zr, = Z - S = Z 



( 10 ) 


One can see, when Zc Z Z*, one gets S Z and 
Zn ^ Z, meaning the charge renormalization is weak. 
However, a,t Z ^ Z*, according to Eqs. ^ and (101, 
we get Zn S Z*. The potential contributed by 
electrons is no longer perturbative. This brings us to the 
new regime of strong renormalization of charge. 


We show that dX Z ^ Z* the net charge Z„e sat¬ 
urates at the level of Z*e. Indeed, when Z grows be¬ 
yond Z*, Zn can not go down and therefore can’t be 
much smaller than Z*. At the same time it can not 
continue going up, otherwise as follows from Eqs. 0 
and ^ with Z replaced by Zn ^ Z*, the total electron 
charge surrounding the charge ZnC at r > R would be¬ 
come Se ~ Zne{Zn/Z*y^'^ ;g> Z„e leading to a negative 
charge seen from infinity. Thus, at Z ^ Z*, the net 
charge saturates at the universal value of the order of 
Z* as is shown in Fig. [TJd. As we emphasized in Fig. [l] 
this result is qualitatively similar to the one obtained for 
heavy nuclei and donor clusters in Weyl semimetals and 
narrow-band gap semiconductors in Ref. 1131 

In the following subsection, we show how the renormal¬ 
ization of charge at Z Z* is realized through certain 
distribution of electrons, in which a structure of “double 
layer” (see Fig. plays an important role. 


C. Radial Distribution of Electrons 

At Z :$> Z*, the charge renormalization is strong and 
the most of the sphere of radius R is completely neutral¬ 
ized by electrons. In the neutral center of the sphere, 
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FIG. 2. (Color online) Potential energy of electrons U{r) = 
—e4>{r) as a function of radius r. <j)Q is defined as n(<)>o) = no, 
where no = 3ZI4ttR^, n(r) is a function of (/>(r) given by Eq. 
Q. The thick solid line (blue) represents the potential profile 
of a cluster of charge Z < Z* which is in the regime of weak 
charge renormalization. The thin solid line (red) represents 
the potential of a cluster at Z ^ Z* in the strong renor¬ 
malization regime, where the two vertical dotted lines show 
edges of the “double-layer” structure of width ~ d <C I?. The 
horizontal dashed line (black) indicates the position of the 
chemical potential p = 0. r* is the external radius of the 
collapsed electron gas where Thomas-Fermi approach fails. 


the electron density n(r) = no, where no = iZ/ATvR^ is 
the density of the positive charge inside the cluster. The 
corresponding “internal” electric potential 0 i„(r) = (j)Q 
where (j)o is given by n(0o) = no using Eq. 0 ■ ^m(^) 
is then oc (noa^)^/^ oc [Z/(i?/a)^]^/^. Outside the clus¬ 
ter, when the charge is renormalized to .Z„, one gets a 
potential tpoutiT) similar to Eq. Q with Z replaced by 
Zn- Since Z^ is ~ Z* where Z* is given by Eq. (|^, 
we get (/)outir) cx at a distance r of the order 

R. Thus, close to the cluster surface, the ratio of the 
outside potential (poutir) to the inside potential 0 i„(r) is 
~ {Rlaf/yZ^/^ < 1 since Z ^ Z* - {R/af/\ This 
indicates a sharp potential drop across the sphere surface. 

At 0 < R — r R, there’s a thin layer of uncom¬ 
pensated positive charges. At0<r — a higher 

potential than farther away means a larger electron con¬ 
centration that forms a negative layer close to the sur¬ 
face. This “double-layer” structure resembles a capacitor 
which quickly brings the potential down across the sur¬ 
face as shown in Fi g. An analogous structure also 
exists in heavy nucleP^^with charge Z ^ 

To make the analysis more quantitative, one needs to 
know the specific potential profile in this region. Using 
Eq. (§, we get the general equation of (()(r) in the spher¬ 
ical coordinate system: 


( 


d 

dr 


+ 



( 


_d 

dr 


+ 



[n{r) - no] , r < i? ( 11 a) 

A^/^e 

-^/^n(r), r > R (lib) 


Near the cluster surface, we can approximately use a 
plane solution of <^(r), i.e., ignore the 2 /r term on the 
left side. This kind of solution for r > i? is already 


knowrfi^ 


Hr) = 


-16/7 
b \a) 


Co e 


X + d 


8/7 


( 12 ) 


where x = r — R<^R is the distance to the surface 
and d R is the characteristic decay length to be de¬ 
termined, C 3 « 6 . Correspondingly, the radial electron 
concentration at r > i? is given by 


n(r)r^ = 


C 4 


1 nx 


^3/7 1)3 j 

C 4 1 / ^ \ 
^ 3/753 I a/ 


24/7 

( 

/ a ^ 

12/7 

1 

yx + d ^ 

I 

24/7 

^ a ^ 

12/7 

1 

yx + d j 

) 


(13) 


R^ 


where r « i?, C 4 « 1 . 

Since the “double-layer” structure resembles a plane 
capacitor, near the surface, the potential drop is nearly 
linear with the radius. Using Eq. (12), one can get 
(f>{r) « (1 — 8x/7d)(j){R) atO <x = r — R<^d, 
which gives the electric field 8/()(i?)/7d inside the “dou¬ 
ble layer”. At r < i?, this electric field persists and gives 
(j){r) «(!-!- 8x/7d)4>{R) at Q < x = R —r <€. d. As r fur¬ 
ther decreases, the positive layer ends and the potential 
crosses over to the constant value (po given by n{(l)o) = no 
using Eq. ([^. This boundary condition gives 


d = 


C 5 


1/4 


A^/^\aJ (noa3)7/i2’ 


(14) 


where C 5 « 2. By expressing no in terms of Z and i?, we 
get d/R oc {Z*IZ)'^'^‘^ < 1 at Z > ZA 

According to Eq. (12), when x is comparable to R 


and the plane approximation is about to lose its validity, 
4>{r) is oc (i?/a)“®/^. It is weak enough to match the 
low electric potential (jjoutir) cx (i?/a)“®/^ caused by the 
renormalized charge ~ Z* at r ~ R. The plane solu¬ 
tion then crosses over to the potential (j)out{r) oc Z*^jr^ 
which is the asymptotic form at large distances. 

A schematic plot of the potential energy U{r) = 
—e4>{r) as a function of radius r is shown in Fig. iby 
the thin solid line (red). The corresponding radial distri¬ 
bution of electrons is shown in Fig. by the thick solid 
line (red). 


So far, we have got a l/r® potential (/)(r) and l/r^^/s 
radial electron concentration n(r)r^ at r i? in both 
weak and strong charge renormalization cases. However, 
as the electron density decreases to certain extent so that 
the Fermi wavelength A is comparable to the radius r, 
the gas is no longer degenerate and the Thomas-Fermi 
approach fails. Since A ~ n{r)~^^^, we get this radius 
r* as Z* a at Z ^ Z*. One should then return to the 
Schrodinger equation used for a single electron. Since the 
uncertainty principle estimates that the kinetic energy 
decays as 1 /r^ while the potential energy is cx —Z*^jr^, 
the potential energy is smaller than the kinetic energy in 
magnitude at r > r*, which means electrons can not stay 
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FIG. 3. (Color online) Radial electron concentration n{r)r^ 
as a function of radius r. The thick solid line (red) represents 
the inner collapsed electrons at r < rc where the dielectric 
response is nonlinear. The thin solid line (blue) shows the 
electrons belonging to the outer shell which form the standard 
Thomas-Fermi atom with the renormalized nucleus of charge 
Z* at r > rc, where the dielectric response is linear. This 
electron gas ends at the Bohr radius ta = while most of 
them are at radius rm = K,b/ . The dashed line (black) 
denotes the electrons forming a Thomas-Fermi atorrP^with a 
nucleus of charge Z when Po is infinity and there’s no range 
with nonlinear dielectric response. The reduction of electron 
density in the outer shell of electrons due to the collapse is 
substantial. The reason this is not immediately seen from the 
difference of height between the dashed line (black) and the 
thin solid line (blue) is that we use a logarithmic scale here, 
no is defined as 3Z/47rP®. This graph is plotted at b — 0.35 A, 
a = 3.9 A, A — 0.9, R — 4.4a, k = 20000, no = 0.8/a®. 


at radii larger than r*. One can also find that using the 
Thomas-Fermi solution (jjautir), the total electron num¬ 
ber calculated at r > r* is ^ 1 , which again indicates 
there’s no electron at r > r* considering the discreteness 
of electron charge. As a result, the tail of radial 

electron concentration will not continue to infinity but 
stop at radius r*. This is a semi-classical result. Quan¬ 
tum mechanical analysis shows that the electron density 
does not go to zero right at r* but decays exponentially 
after this point. Since this decay is fast and brings very 
small corrections to the end of the inner electron gas, we 
don’t consider it here. 

At K = oo, the rest of the electrons are at the infinity 
so that we are dealing with a positive ion with charge Z*. 
At finite but very large k, at certain distance from the 
cluster, the field is so small that P > Att /tiAPo is no 
longer satisfied and the linear dielectric response is recov¬ 
ered. Things then become quite familiar. Electrons are 
mainly located between ri = KhjZ* and ta = i^b with 
the majority at radius = nb/Z*^^^ as given by the 
Thomas-Fermi modeP^. Although quantum mechanics 
gives a nonzero electron density at r < ri, the num¬ 
ber of total electrons within this radius is only ^ 1 and 
can be ignored. So approximately, when ri 3> r*, i.e., 
K » (Z*)^ ~ there’s a spatial separation be¬ 

tween inner collapsed electrons and outer ones that form 
the usual Thomas-Fermi atom with the renormalized nu¬ 
cleus. When K is not so big, such separation is absent. 


which actually happens more often in real sitnations. The 
inner tail then connects to the outer electrons with the 
Thomas-Fermi approach valid all the way and the dielec¬ 
tric response becomes linear at r = Cc oc . 

One should note, as long as k is large enough to satisfy 
r^n ^ R which gives k ^ (ii/a)^°/’^, the majority of 
the outer electrons located at haven’t intruded into 
the cluster or the highly screening double-layer structure 
near the cluster surface. The charge renormalization pro¬ 
cess remains undisturbed and the total net charge seen 
by outer electrons is still Z*. The corresponding radial 
electron concentration n(r)r'^ is shown in Fig. At 
K < (i?/a)lo/^ in most of the space the dielectric re¬ 
sponse is linear. In that case, almost all electrons re¬ 
side in the cluster with only some spill-over near the sur¬ 
face. The positive and negative charges are uniformly 
distributed inside the cluster as described by the Thomp¬ 
son “jelly” model. 


III. CYLINDRICAL DONOR CLUSTERS 


In some cases, the donor clusters are more like long 
cylinders than spheres. For this situation, a cluster is 
described by the linear charge density rje while its radius 
is still denoted as R. We use a cylindrical coordinate 
system with the z axis along the axis of the cylinder 
cluster and r as the distance from the axis. We show that 
when the charge density rje is larger than certain value 
rjcC, electrons begin to collapse into the cluster and the 
charge density is weakly renormalized. When rj exceeds 
another value rj* 3 > 77c, the renormalization becomes so 
strong that the net density remains ~ rj* regardless 
of the original density 77 (see Fig. [^. Our problem is 
similar to that of the charged vacuum condensate near 
superconducting cosmic string^^, and is also reminiscent 
of the Onsager-Manning condensatiorP^^ in salty watei!^. 

Renormalization of Linear Charge Density. — For a 
uniformly charged cylindrical cluster with a linear charge 
density rye, similar to what we did in Sec. II, we get 
D{r) = 2 ri{r)e/r, where r]{r) = rjr'^/R'^ at r < R and 
r]{r) = T] at r > R. We then can calculate the electric 
field using Eqs. ([^ and ([^ and get the electric potential 
(j){r) as: 


cj){r) 

cj){r) 


A zryeN^/S 1 

1^\^) [aW 

A r rje 1 


AR^j 


, 0 < r < i? 


R < r 


(15a) 

(15b) 


with (j)[r = 00) chosen to be 0. The corresponding po¬ 
tential energy U{r) = —e(j){r) is shown in Fig. by the 
thick solid line (blue). Using the Schrodinger equation 
and setting the momentum p ~ H/ 2 r, we find that the 
tightly bound states of electrons, in which electrons are 
strongly confined within the cluster (at r < i?), exist only 
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when r] > T]c where 


'i) 


1/3 


(16) 


which, contrary to the Zc value got in the spherical case, 
does not depend on R. Electrons begin to collapse into 
the cluster at ly > t/c and in the beginning they are lo¬ 
cated near the axis where the potential energy is lowest 
(see Fig. [^. With increasing ry, the electron density 
grows and one can adopt the Thomas-Fermi description. 
Using Eq. 0 and (|15b[), one gets the electron density 


(r) oc 1/r^ at r > i? and the total number of collapsed 
electron per unit length is 

0 = y n{r)2'Krdr = 0.5ry ^ (17) 


where 


ry = 




1/3' 


9/7 


a) 


(18) 


The net charge density ?y„e is then renormalized to 


r]n = ri-e = r] 


1 - 0.5 



(19) 


At ly <?; ry*, the renormalization of charge density is 
weak and ry„ grows with ry. At ry > ry*, the number of col¬ 
lapsed electrons is large and the renormalization effect is 
strong. Most of the cluster is then neutralized by elec¬ 
trons and the final net density ry„ is much smaller than ry. 
Following the logics similar to those in the spherical case, 
and by using Eq. 0 , one can show that ry„ reaches a 
saturation value of ry* at ry 3> ry*. The dependence of ry„ 
on ry is shown in Fig. which resembles Fig. 



ri’iR) 


FIG. 4. (Color online) Number of collapsed electrons per unit 
length 9 and renormalized net charge density ry„ as a function 
of cluster charge density ry. The thick solid line (blue) shows 
ry„(ry). The thin solid line (red) represents ^(ry). The dashed 
line (black) is a guide-to-eye with ry„ = ry. 6(r]) oc ry®A at 
rye <C ry < ry*. 

Radial Distribution of Electrons. — At ry ry*, there 
are lots of collapsed electrons inside the cluster where 


n(r) = no = ry/rri?^ and the potential energy is low. 
Again, there’s a “double-layer” structure on the surface 
that provides steep growth of potential energy with r at 
r = R. Close to the cylinder surface at 0 < r — i? ^ i?, as 
for the sphere, we can approximately use a plane solution 
of (f>{r) as given by Eq. (12). The expression of the 
characteristic decay length d is also the same as in Eq. 
(14). When a; = r — i? is comparable to i?, the plane 
solution crosses over to the fast decaying potential oc 1 /r^ 
as given by Eq. ( |15b ) with ry replaced by ry„ ~ ry*. A 
schematic plot of the potential energy U{r) = —e(j){r) is 
shown in Fig. 



FIG. 5. (Color online) Potential energy of electrons U{r) = 
—efir) as a function of radius r. (j>Q is defined as n{(j)Q) = 
where no = r)/-KR^, n{r) is a function of ^{r) given by Eq. ( 1 ^. 
The thick solid line (blue) represents the potential profile of a 
cluster of charge density ry < ry* which is in the regime of weak 
renormalization of charge. The thin solid line (red) represents 
the potential of a cluster with ry ry* which is in the strong 
renormalization regime. The two vertical dotted lines show 
edges of the “double-layer” structure of width ~ d <C A. The 
horizontal dashed line (black) indicates the position of the 
chemical potential fi — 0. 

This potential produces a universal tail of electron den¬ 
sity n(r) ^ 1/r^. The corresponding radial electron con¬ 
centration n{r)r is ^ 1/r^. Since the Fermi wavelength 
A ~ n(r)“^/^, we get A ~ r, i.e., the Thomas-Fermi ap¬ 
proach is only marginally valid. The collapsed electrons 
extend until the linear dielectric response is recovered 
and then connect to the outer electrons. 


IV. FINITE-TEMPERATURE IONIZATION OF 

CLUSTER ATOMS AND ITS EXPERIMENTAL 
IMPLICATIONS. 

So far, we dealt with very low temperature. At a finite 
temperature T, the neutral cluster atom can get ionized 
due to the entropy gain of ionized electrons. The donor 
cluster atom becomes a positive ion with charge Zi{T)e. 
Our goal below is to find this charge. 

We assume that we have a small but finite three- 
dimensional concentration N of spherical clusters and 
the charge Zi{T) < Z* , i. e., the outer electron shell 
is still incompletely ionized. Such a cluster can bind 
electrons with an ionization energy Zi{T)'^e^/K^b. We 
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can find Zi(T) by equating this energy with the decrease 
in the free energy per electron kBTla{nQ/n) due to the 
entropy increase, where ks is the Boltzmann constant, 
n = Zi{T)N is the concentration of ionized electrons 
and no = 2/A^ with A = ^2TTh?/m*kBT as the De- 
Broglie wavelength of free electrons at temperature T. 
At K = 20000, b = 0.29 A, m* = 1.8me where rUe is 
the electron mass and N = 10^® cm“^ (estimated from 
that the concentration of total donor electrons is around 
10^® cm“® and each cluster contributes ^ 300 donor elec¬ 
trons), we get Zi{T) > Z* a.t T > 8 K with Z* = 100 
which is a reasonable estimate. This shows that the 
outer electrons are completely ionized at temperatures 
that are not too low. For the inner core electrons, the 
dielectric response is nonlinear and the attractive po¬ 
tential is stronger. So the ionization energy is higher 
~ A{Z* e/A tt)^ / 5Pq for electrons air ~ R. At i? = 4a, 

it is found that only at T > 450 K can the inner elec¬ 
trons be ionized by a considerable quantity (the 
tail is completely stripped then). So the inner electrons 
are robust against the thermal ionization. 

For the cylindrical cluster, since it can effectively be 
regarded as the assembly of sphere clusters, one can ex¬ 
pect that electrons are harder to be stripped off than in 
the spherical case. The thermal ionization is thus some¬ 
what weaker. For the outer electrons, the dielectric re¬ 
sponse is linear and the ionization degree is determined 
by the Onsager-Manning linear density 770 It can 
be derived as follows. The potential energy of electrons 
caused by the cylindrical charge source with a linear den¬ 
sity rjie grows with the radius r as 2rjie^ \a[r/r ^)/k while 
the entropy increases as 2k b ln(r/ro), where rp is a cho¬ 
sen reference point. In equilibrium, by equating the en¬ 
ergy increase and the entropic decrease of the free energy 
one gets the critical concentration 77 ^ = 770 M = nkBT/e^. 
This is the universal value of the net charge density, 
which depends only on the temperature T and the di¬ 
electric constant k of the media similar to the case of 
DNA in salty wateJ^. 

When the outer electrons are completely ionized, the 
charge density is expected to be T]*e given by Eq. (18). 
Taking R = 4a, A = 0.9, k = 20000, a = 3.9 A and 
b = 0.29 A, we get r]oM ^ 77 * at T > 10 K. Thus, the 
ionization of the outer shell proceeds until T grows to 
10 K. Since the product kT is almost fixed at T > 10 K, 
VOM actually stops growing with the temperature at T > 
10 K. This practically means that only outer electrons 
are ionized at T > lOK. The inner electrons are well 
preserved against the ionization even at T > 10 K. 

Thus, in both spherical and cylindrical cases, the outer 
electrons are thermally ionized at not very low temper¬ 
ature while the inner ones are mainly kept by the clus¬ 
ter. The final observable charge or charge density is then 
equal to Z* or 77 *. Below we will discuss how one can 
observe the thermal ionization. 


Experimentally, charged clusters can be created con- 
trollably on the surface of LAO/STO structure when the 
LAO layer is of subcritical thickness < 3 unit celP. A 


conducting atomic force microscope (AFM) tip is placed 
in contact with the top LaAlOa (LAO) surface and biased 
at certain voltage with respect to the interface, which is 
held at electric ground. When the voltage is positive, a 
locally metallic interface is produced between LAO and 
STO where some positive charges are accumulated in the 
shape of a disc. The same writing process can also create 
a periodic array of charged discs. 

Let’s first concentrate on a disc of positive charge cre¬ 
ated in this manner on the STO surface. Close to the 
surface and in the bulk STO, one should apply the plane 
solution given by Ref. [18] and repeated by Eq. (12) 
above. When the distance r from the disc center is large, 
i.e., r ^ R, the disc behaves like a charged sphere. Our 
results for a sphere are still qualitatively correct in this 
case. 

In a periodic array of highly charged discs with period 
2 L, the linear concentration of free electrons responsi¬ 
ble for the conductance at a very low temperature is of 
the order of n(L)L^, where n(r) is the electron density 
around a spherical donor cluster given by Sec. 11. When 
the overlapping parts between neighboring discs belong 
to the outer electron shells, the corresponding density at 
r = Lis that of a Thomas-Fermi atom with charge Z*. In 
this situation, the overlapping external atmosphere forms 
conductive “bridges” between discs at low temperature. 
When T increases, however, the outer electrons are ion¬ 
ized and the bridges are gone. These free electrons spread 
out over the bulk STO. At T < 30 K, electrons are scat¬ 
tered mainly by the Coulomb potential of donors and 
the corresponding mobility decreases with a decreased 
electron velocity. For the electrons ionized into the vast 
region of the bulk STO, they’re no longer degenerate, 
so their velocity becomes much smaller at relatively low 
temperature. This results in a much smaller mobility of 
the ionized electrons than those bound along the chain. 
Their contribution to the conductivity is thus negligible. 
The system becomes more resistive due to the ionization 
and one can observe a sharp decrease of the conductivity 
along the chain. It may be possible for the chain to tran¬ 
sition from a conducting “line” to multi-quantum-dots. 
Interesting phenomena of the conductance behavior such 
as the Coulomb blockade can probably emerge. 


V. CONCLUSION 


In this paper we have studied the structure of a many- 
electron “atom” whose center is a strongly charged donor 
cluster. It is determined by the collapse of electrons to 
the cluster in SrTiOs due to the nonlinear dielectric re¬ 
sponse at small distances from the cluster surface. For a 
spherical cluster, when its charge number exceeds a crit¬ 
ical value Zc, the potential well inside it becomes deep 
enough to trap electrons despite their high kinetic en¬ 
ergy. In the beginning, the cluster charge Ze is weakly 
renormalized by the electrons. When Ze grows beyond 
another value Z*, the number of collapsed electrons be- 





comes large so that the most of the cluster is neutral¬ 
ized and the charge is strongly renormalized to Z*. This 
strong renormalization is realized via a “double-layer” 
structure on the cluster surface. The corresponding po¬ 
tential profiles and radial distributions of electrons are 
investigated. The critical and saturation values are found 
to be dependent only on the cluster radius: Z^ ~ (^/o) 
and Z* ~ At zero temperature, a renormal¬ 

ized cluster with charge Z* is the nucleus of a Thomas- 
Fermi atom. This nucleus is surrounded by the exter¬ 
nal electron atmosphere which is sparse due to the weak 
Coulomb interaction at a large dielectric constant. These 
outer electrons can easily be stripped off by the ther¬ 
mal ionization leaving only the compact ionic core with 
charge Z*. The case of a cylindrical donor cluster is dis¬ 
cussed as well where similar results are found. Namely, 
when its linear charge density T^e is larger than r\ce ~ e/a, 
electrons start collapsing to the cluster. At 77 3> ry* where 


is the saturation density, the net charge 
density of the cluster ry„e stays at the level of rfe. At 
zero temperature, this renormalized cylindrical cluster is 
surrounded by a sparse electron atmosphere which can be 
ionized at temperatures larger than 10 K. We also dis¬ 
cuss how one can verify our predictions by measuring the 
conductivity of a chain of charged discs on the surface of 
LAO/STO structures. 
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